Theoretical model for the structural phase transition at the metal-insulator transition 

in polymerized KCgo 
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The recently discovered structural transition in polymerized KCgo at about 50 K results in a 
doubling of the unit cell volume and accompanies the metal-insulator transition. Here we show that 
the (a + c,b, a — c) superstructure results from small orientational charge density waves along the 
polymer chains and concomitant displacements of the surrounding K + ions. The effect is specific 
for the space group Pmnn of KCeo and is absent in Rb- and CsCeo (space group 72/m). The 
mechanism is relevant for the metal-insulator transition. 
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I. INTRODUCTION 



Among the alkali-metal-doped fullerides A^Cgo, where 
A=K, Rb, Cs, the x = 1 compounds Q have at- 
tracted attention, because they form plastic cubic crys- 
talline phases with rock salt structure at high temper- 
ature (T > 350 K) and stable polymeric phases g of 
reduced symmetry at lower T. In addition, via rapid 
quenching, a metastable dimer phase was obtained §|, |. 
It was suggested that the dimer structure should be a 
Peierls insulator ||| . In the following we will restrict our- 
selves to the polymer phases. There, the Ceo molecules 
are linked through a [2+2] cycloaddition and form chains 
along the former cubic [110] direction. The orientation of 
the polymer chains is characterized by the angle vb of the 
planes of cycloaddition with the former cubic [001] direc- 
tion. The structure of KCeo is orthorhombic || |(| [space 
group Pmnn, Fig. 0(a)] ; the orthorhombic a axis is paral- 
lel to the axis of polymerization and the c axis parallel to 
the former cubic [001] direction. The Pmnn structure is 
characterized by alternating orientations ±ip of the poly- 
mer chains in successive (a, b) planes and the same ori- 
entation of the chains within one plane, « 45°. On 
the other hand, the structure of both polymerized RbCgo 

fand CsCgo @ is monoclinic, space group 12 /m [Fig. 
b)]. Here the polymer chains have the same orienta- 
tion tp not only within the same (a, b) plane but also in 
successive (001) planes. 

In addition to the structural differences, the electronic 
and magnetic properties of KCeo on one hand and Rb-, 
CsCeo, on the other hand, are different (for a review, see 
Rcf. ||). Neither the differences in structure nor those 
in electromagnetic properties can be simply related to 
the relatively small differences in the lattice constants 
among the ACgo compounds. We then conclude that 
more alkali-metal-specific effects are relevant. To explain 
the structural differences of the polymer phases of KCeo 
(space group Pmnn) and Rb- and CsCeo (space group 
72/m), it has turned out that the quadrupolar polar- 
izability of the alkali-metal ions is the decisive alkali- 
specific-characteristic |qJ, 
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ability of the alkali-metal ion A + is related to the aver 



FIG. 1: Crystal structures projected onto the crystallo- 
graphic (b,c) plane: (a) Pmnn, (b) 72/m. The thick bars 
represent the projection of the cycloaddition planes. Poly- 
merization occurs along a. The alkali-metals located in (b, c) 
planes and at ±o/2 are denoted by full and empty circles. 



age radius dx of its first valence electron d shell. Indeed, 
d-Rh = 1-82 A and des = 1-87 A are close to each other 
but quite different from c?k = 1-47 A. The partial occu- 
pancy of the first excited d state in the solids is possible 
because of the large interstitial space available for the 
alkali-metals, this relatively large space being a unique 
property of the ACgo compounds. 

Recently, a combined study of the electronic and struc- 
tural properties of KCeo has been carried out by Coulon 
et al. pdj ] : x-ray diffraction studies have revealed a struc- 
tural phase transition in polymerized KCeo at 60 K 
< T c < 65 K, while ESR measurements have shown the 
existence of a metal- insulator transition at T w 50 K, 
which stresses once again the close connection between 
structural and electronic properties in the ACeo alkali- 
metal fullerides. Above the structural critical temper- 
ature T c , the Bravais lattice of polymerized KCeo is or- 
thorhombic; at T c , the crystal changes to a (o + c, b, a — c) 
superstructure, which can be viewed as a doubling of the 
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basis vectors a and c. Concerning the possible nature of 
this structural phase transition (i.e., the doubling mech- 
anism), it was pointed out by Coulon et al. [[fl| (i) that 
displacements of the Ceo centers of mass can be ruled out 
and (ii) that a combination of a charge modulation on the 
Ceo monomers and large correlated K displacements is an 
appealing hypothesis. 

We note that such a combined structural and metal- 
insulator phase transition has not yet been observed in 
RbCeo [fl2|. Bearing in mind the structural differences 
of the polymer phases of KCeo and RbCeo, the pres- 
ence/absence of a second structural phase transition in 
AC 6 o seems to be related to the precise orientations of 
the polymer chains in Pmnn and 12 /m, respectively. 

Here we propose a doubling mechanism that accounts 
for the observed structural phase transition in KCeo on 
one hand, and for the absence of such a transition in 
RbCeo on the other hand, both in agreement with the 
present experimental knowledge of these compounds. 
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II. DOUBLING MECHANISM 

We consider rigid KCeo and RbCeo crystals in the poly- 
merized phase. For KCeo, this phase has space group 
Pmnn, while for RbCeo the space group is 12 /m. The 
12 1 m space group has a monoclinic Bravais lattice; we 
will treat it as an orthorhombic lattice, however, since the 
"monoclinic" angle a — 90.316° of RbCeo obtained by 
neutron scattering measurements [|l3| is extremely close 
to 90° . The centers of mass of the N Gqq monomers are 
then located on the lattice points of a body-centered or- 
thorhombic lattice with basis vectors a = aex, b = bey, 
and c = cez, where (ex,ey,ez) are the basis vectors of 
the underlying Cartesian coordinate system. The axes of 
polymerization are parallel to a. A lattice point is labeled 
by indices ft = (ni, ri2, TI3), which are either all integer 
numbers or all integer numbers +5, corresponding, re- 
spectively, to the corner and the center points of the unit 
cells. The position vector of lattice point n reads 



X{n) = n\d + n^b + n 3 c. 



To each polymer chain, a rotation angle ip can be as- 
signed. We let ip = correspond to the situation where 
the polymer chain is in the standard orientation, which is 
defined as the orientation where the plane of cycloaddi- 
tion is parallel to the (a, c) plane. The angle ip then mea- 
sures a counterclockwise rotation of the polymer chain 
about a. More generally, the orientation angle can be 
seen as a property of a single Ceo monomer; therefore we 
write ip = ip{n). However, since all the Ceo monomers in 
the same polymer chain have the same orientation angle 
and since a polymer chain can be addressed by the indices 
(712,713), ip is independent of the index n\. Furthermore, 
the Pmnn structure is characterized by an alternation of 
the orientations of the polymer chains along the c axis 
only; hence for KCeo ip is also independent of the index 



ip{n) = ip(n 3 ) 
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^KCso- 



(2.2) 



The structure 12 /m is characterized by an equal orienta- 
tion of all the polymer chains in the crystal. Therefore, 
one has for RbCeo: 



ip(n) = ip 



RbO 



(2.3) 



The chain orientation angles have been redetermined re- 
cently by neutron scattering experiments jl3fl : t/'kCbo = 
50°, ^ Rb c 60 = 46°. 

Due to the charge transfer of one electron from an 
alkali-metal atom to a Ceo molecule, a charge distribu- 
tion p with total charge — e exists on every Ceo monomer. 
For a given Ceo monomer n, we introduce a Cartesian 
coordinate system with basis vectors {e x ,e y ,e z ) and the 
center of mass of the monomer as origin. The charge 
distribution can then be written as 



P = p(r; ip{n)) = p(x, y, z; tp(n)), 



(2.4) 



where r = xe x + ye y + ze z . 

As a mechanism for the doubling of the lattice basis 
vectors a and c, we suggest the following: while retaining 
the rigid structure for the C nuclei and the closed tt- and 
a -shell electrons, small orientational deviations of the va- 
lence electron density p on every Ceo monomer from this 
structure are allowed to occur, in such a way that these 
deviations (and therefore equivalent lattice points) have 
periodicities 2a, b, and 2c. A rotation of the electron 
density can occur because there are orbital degrees of 
freedom for one valence electron on the threefold degen- 
erate t\ u lowest unoccupied molecular orbital (LUMO) 
level [Q. The valence electron density depends on the 
coefficients of expansion in terms of these three orbitals 
and this leads to its effective rotation. We will limit our- 
selves to rotations of the charge distributions about the 
polymer chain axes. To denote the angular deviation of 
the charge distribution from ip(n), we use the notation 
Aip(n): 

p = p(f; ip(n)+Aip(n)) = p(x, y, z; tp(n)+Atp(n)). (2.5) 
(2-1) In order to have the desired periodicity changes, we take 



Atp(n 1 ,n 2 ,n 3 ) = Aip(m + ^, n 2 + ^, n 3 + ^) 
+A-0O if n\ +773 even 



— A-0o if 77i + 713 odd, 



(2.6) 



with 77i,772,7i3 e Z. In Eq. (2.6), Ai/>o represents the 
angle measuring the deviation from ip. The mechanism 
is illustrated schematically in Fig. ^. In Sec. IV we will 
comment on the charge distribution model to be used. 

It is convenient to introduce spherical coordinates 
(r, 9, (p), related to the Cartesian coordinates (x, y, z) via 



y = r sin a cos < 
z = r sin 6 sin <; 
x = rcos9, 



(2.7) 
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FIG. 2: Illustration of the proposed doubling mechanism. Shown is a projection of the Ceo molecules, symbolized by circles, 
onto the (b, c) crystallographic plane (compare with Fig. [j]) . The alkali-metal atoms have been omitted for clarity. The solid 
and dashed lines in the circles represent the orientations of the charge distributions: (a) KCeo, (b) RbCeo- The solid lines 
correspond to the original Pmnn (a) and 12 /m (b) structures, the suggested doubling mechanism consists of angular deviations 
of the electronic density distributions as indicated by the dashed lines. The proposed scheme results in a doubling of the lattice 
vector c, while the lattice vector b remains the same. The original (projection of the) unit cell is shown as a solid frame, the 
(projection of the) unit cell taking into account the electronic deviations is shown as a dashed frame. Along the a axis, a 
doubling pattern (angular deviations of the charge distributions) similar to the doubling along the c axis occurs. 



since a counterclockwise rotation of the charge distribu- 
tion p over an angle ip + Aip about the a axis is then 
simply achieved by replacing <fi by <fi — ip — Aip: 



p{r, 9, 4>; ij){n) + Ai/>(n)) = p{r, 6,<f>- i>(n) - A^(n); 0). 
We simplify the notation and write from now on: 



III. POTENTIAL ENERGY 



In order to determine whether the above-described 
electronic density distortions can occur or not, it is nec- 
essary to investigate the potential energy of the crystal. 
We consider the contribution U to the potential energy 
arising from Ceo-Cgo interactions, which can be written 
as a sum of pair potentials: 



p(r,M) = p(r,M;0). 



(2.9) 



U 



(3.1) 
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Since we are interested in how the potential energy of the 
crystal is influenced by the above-described orientational 
deviations Aip(ft) of the charge distributions on all Cg 
monomers, only electronic interactions have to be taken 
into account (van der Waals-type interactions refer to 
the neutral C cores, which in our model are assumed to 
constitute a rigid structure). The electrostatic energy of 
two Cg monomers at lattice points ft and ft' is given by 



U(n,n') 



1 



4"7re 



dr I dr 



p(r;a(ri))p(r';a(n')) 
\r — f — X(ft' 



(3.2a) 



where we have introduced the notation 



a(ft) = ij)(n) + Aip(ft). 



(3.2b) 



The integration variables r and r* in Eq. (3.2a) refer to 
the local coordinate systems associated with the respec- 
tive lattice sites ft and ft'; hence the appearance of the 
relative position vector X(n' — n) — X(n') — X(n). Using 
the previously introduced spherical coordinates, U(ft,ft') 
can be rewritten as 



U(n,n) = 



1 



dr / sin 



oo p 7T />2n 

r a dr' \ smO'de' 



47T60 Jo JO JO JO JO J" 

p(r, 6,(f>- ip{ft) - Ai){n))p(r',0', tj>' - tp(n') - A^(n')) 



\f — f" — X(ft' — ft)\ 

\f — f 1 — X(ft' — ft)\ = { [r cos 9 — r' cos 9' — {n'-^ — ni)a] 2 + [r sin 9 cos 4> — r' sin 9' cos 4> — (n' 2 — ri2)6] 2 

1 /2 

+ [r sin 9 sin <f> — r' sin 9' sin <p' — (n' 3 — ?i3)c] 2 } 

I 



(3.3a) 



(3.3b) 



As will be discussed in more detail in Sec. IV, the 



charge distribution p can be expanded in even multipolcs 
I = 0, 2, 4, . . . . The lowest-order term containing angular 
dependence will therefore be a monopole-quadrupole in- 
teraction. As a consequence, the angular dependent part 
of the potential energy U, being of the Coulomb-type, 
decreases fast enough with the distance, and it is jus- 
tified to consider a limited number of nearest-neighbor 
interactions. We write Eq. (pjj) as 



U 



where 



V(rt) = J2u(rt.,ri + p). 



(3.5) 



The summation over p runs over the fourteen nearest- 
neighbor sites with relative indices p — (±1,0,0), 
(0,±1,0), (0,0, ±1), (±i,±i,±i). An analysis of the 
occuring orientations of the charge distributions on the 
various Cg~ monomers reveals that only two different 
types of lattice sites exist: if we let m, n-z, and nz 
be integers, then all sites with indices [n\, 712,713) or 
(ni + ^ } n.2 + 5,^3 + 5)) satisfying ri\ + even, have 
an equivalent neighborhood. We call these sites type I 
sites. The other sites {n\ + n 3 odd), which we label as 
type II sites, have also an equivalent environment, but 
one that is different from the type I neighborhood. We 



summarize this observation by writing 

F(ni, 712,713) = V(n\ + i,n 2 + 7^3 

V 1 if m + 713 even 
V u if ni + 113 odd, 



(3.6) 



with ni,7i2,7i3 G Z. Note that Eq. (3.6 ) is consistent 
with the imposed periodicity conditions (|2.6|) . The po- 
tential energy of the entire crystal due to all electrostatic 
CgQ-CgQ interactions is obtained by carrying out the sum- 



(3-4) mation of Eq. (3.4), which runs over all lattice points: 



U = ^(V* + VK). 



(3.7) 



The functions V 1 a nd V 11 depend only on the angle A-0o , 
introduced in Eq. ( |2.6| ), which is a measure for the devi- 
ation from the undistorted structure. To emphasize this 
dependence, we write 



U=U(Ai/j q ) 



IV. CHARGE DISTRIBUTION 



(3.8) 



The key point in the proposed model is the expression 
for the electronic density p{r) of a Cg~ unit in a poly- 
mer chain. Polymerization leads to a reduction of the 
symmetry (in comparison with the icosahedral symme- 
try of Cgo) and the charge distribution can be expanded 
in multipoles with I = 0, 2, 4, . . . . In earlier theoretical 
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work [T^] , a quadrupolar charge distribution model was 
used in explaining the structural phase transition from 
the cubic (unpolymerized) to the orthorhombic (poly- 
merized) phase of the ACeo alkali-metal fullerides. By 
using a simplified Slater-Koster tight-binding approach 
to determine p(r), it has been shown |H| that such a 
quadrupolar model is a reasonable first approximation. 
(Due to the F>2h symmetry of a polymer chain, only even 
multipoles occur.) However, by experimenting with vari- 
ous quadrupolarlike electronic densities, we find that the 
energy U depends very sensitively on the precise location 
of the point charges used in constructing the charge dis- 
tribution and that here no conclusion can be drawn based 
upon quadrupolar models. In our view, it is necessary to 
go beyond the quadrupolar contribution and to take into 
account higher multipoles. Indeed, if one examines the 
expansion of the angular part p a (9, 4 1 ) of the charge dis- 
tribution in terms of spherical harmonics, calculated in 
Rcf. [l5[ one sees the relevance of the higher-order terms. 
In particular, the (/ — 8, m = 8) term seems to play an 
important role. The importance of relatively-high-order 
multipoles is not surprising: for example, orientational 
ordering in solid Cgo (fullerite) is described using molecu- 
lar and site-symmetry-adapted functions (SAFs) belong- 
ing to the manifolds I = 6, 10, 12 pi]. Furthermore, in 
expanding the van der Waals interactions between poly- 
mers in ACeo in terms of SAFs, we have recently discov- 
ered a similar pattern: the term with (I = 8, m = 8) is 
remarkably dominant. 

These observations lead us to use the electronic density 
of Ref. Fill which can be written as 



If a more refined model for p r is used, the main effect 
is a renormalization of the intersite interactions. On the 
other hand, the angular part p a is anisotropic and results 
in multipolc interactions of different sign and magnitude. 
It is given by 



-O.7O699y 5 1,a (0,0) 
-0.637311^(0,0)] 



0.30659K 



3,s, 



(4.1c) 



In Eq. (4.1c), p a is written as the modulus squared of the 
wave function ipa(9, cj>), which is one of the three degen- 
erate ti u functions of a C^~ ion |18| ]. The real functions 
Y, ,s (0, 4>) are the sine spherical harmonics defined in 



Rcf. |19l. Note that Eq. (4.1c) is exact within the frame- 
work of the tight-binding approach. In Fig. ^, the angular 
electron density p a for a CgQ ion in the standard orienta- 
tion is shown, projected onto the (b, c) and (a, c) planes. 
The charge is mainly concentrated in the equatorial re- 
gion (x = or, equivalently, 9 = tt/2), in agreement 
with recent NMR results [g(J ■ As can be seen clearly in 
Fig. ||(a), the four absolute maxima of the charge dis- 
tribution p a coincide with the centers of pentagons [fl2"| , 
which is very reasonable. Indeed, it is known that in 
the neutral molecule, the centers of pentagons are the 
electron-poor regions and any additional negative charge 
experiences minimal repulsion at these centers. A direct 
consequence of this fact is the Pa3 phase of solid Ceo: 
here, an electron-rich region of one molecule (double C- 
C bond) faces a pentagon of its neighbor El]] . 



p{r, 9, 4>) = -ep r {r)p a {9, </>), 



(4.1a) 



where the charge of the electron — e has been factored out 
in order to ensure that the angular part p a is dimension- 
less. The radial part p r is not relevant for our purposes; 
we locate the charge on a sphere with radius R = 3.55 
A, which is the radius of a Cqq molecule [[l7| : 



Pr(r) 



5(r - R) 
~ R 2 ' 



(4.1b) 



V. POTENTIAL ENERGY CALCULATIONS 



The potential energy U of two interacting C 60 
monomers c an n ow be w r itten down expli citly. Com- 
bining Eqs. ( |3.3a| ), ( |3.3b| ), (|4 laj) , and (|4.1b| ), one gets 



U(n,n') = F I sinddO I dcf> / sm9'd6' I d<p' 



jo 



Pa {9, <\> - VKg) - M(n))p a {9', 4/ - VKgQ - A^(n')) 

D(e,4>,9',(t>';n,n') ' { ' 



with 



e 2 



F = — = 167100 K A, (5.1b) 
47re 

D(9,(j>,9',<f>';n,n) = {[Rcos9 - Rcos0' - (n[ - m)a] 2 + [R sin 9 cos (j> - R sin 9' cos </>' - (n' 2 - n 2 )b} 2 

+ [i?sin6»sin0 - R sin 6' sin <j>' - (n' 3 - n 3 )c} 2 } 1/2 . (5.1c) 

I 

The dependence on Aip of the quantity V 1 + V 11 , to (Pmnn structure) and RbCgo (12 /m structure). The lat- 
which U is proportional, is shown in Fig. W for both KCgo 



6 



a 




N 




. 066 

060 
045 
0.40 
0.3S 
030 
0.£S 
0.20 
0.1S 
010 
0.05 




FIG. 3: Angular part p a (8, </>) °f the charge distribution of a C 60 monomer in the standard orientation, plotted as a contourplot 
on a sphere with radius R — 3.55 A, together with the Ceo cage, (a) Projection of p a onto the (a, c) plane, (b) Projection 
of p a onto the (b,c) plane. The coordinates (x,y,z) and (6,4>) are related to each other via Eq. (2.7) and the constraint 
\J x 2 + y 2 + z 2 = R. On both projections, one can see clearly that the charge is mainly concentrated in the equatorial region 
(x — 0). Note the local maxima in the charge density near the C atoms participating in the cycloaddition bonds, and the four 
absolute maxima located at the centers of pentagons (electron-poor regions of the neutral Ceo molecule). 



TABLE I: Lattice constants of K- and RbCeo, taken from 
Ref. (T = 6 K values for KC 60 and T = 5 K values for 
RbCeo), and of CsC 6 o, taken from Ref. |l| (T = 20 K values); 
units of A. 





a 


b 


c 


KC 60 


9.1185 


9.9010 


14.3467 


RbC 60 


9.0887 


10.0843 


14.1583 


CsC 60 


9.0968 


10.1895 


14.1351 



tice constants used in the calculations are given in Ta- 
ble Q. The potential energy for KCgo exhibits a double 
minimum, which implies that a deviation of the charge 
distributions from the Pmnn structure as described by 
the doubling mechanism of Sec. |l| [Fig. ||(a)] is favored. 
The optimal configuration corresponds to a deviation an- 
gle of Aipo ss 13°. On the other hand, for RbCgo, such a 
deviation [Fig. ||(b)] would never lead to an energetically 
more favorable structure since the energy minimum lies 
at Aifj = 0°. 

For completeness, the result for CsCeo is also shown 
in Fig. ^. Since the 12 /m space group angle is the same 
for both Rb- and CsC 60 , t/jRhC eo = V'CsCao = 46° @, H), 
the only difference between RbCeo and CsCgo arises from 
the different lattice constants (Table |). Intuitively, one 
would therefore expect no qualitative difference between 
Rb- and CsCeo- Indeed, Fig. [| confirms the similarity 
between Rb- and CsCgo- Hence our model implies the 
absence of a (a + c,b,a — c) superstructure (and a con- 
comitant metal-insulator transition) for CsCgo- 



At this point, we note that the essential dependence of 
the potential energy IA due to C^q-C^q interactions on the 
distortion angle A^o is introduced in Eq. (2.6). There, 
the deviation angle ip(n) at lattice site n is defined in such 
a way that the experimentally found doubling scheme is 
automatically recovered. However, by assigning to each 
lattice site ft an order parameter S(n) = ±1, correspond- 
ing to a deviation angle A^(n) = ±Aipo, one can show 
rigorously that this doubling scheme has the lowest po- 
tential energy. This analysis is carried out in detail in 
Appendix |Aj. There it is shown that the average order 
parameter (S(n)} is given by 



(S(n))=r,(-l) n ^, 



(5.2) 



where rj is the orde r parameter amplitude. The conden- 
sation scheme (5.2) is indeed equivalent to Eq. ( |2.6| ), ex- 
pressing periodicity doubling along a and c and no change 
in periodicity along b. 



VI. DISPLACEMENTS OF THE 
ALKALI-METAL IONS 

It is clear from Fig. || that the valence electronic density 
deviations Atf>(n) from the Pmnn structure are already 
sufficient to account for the periodicity doubling along 
the a and c axes in KCgo observed by the x-ray scattering 
measurements described in Ref. [ll]. In this section wc 
investigate the role of the alkali-metal ions. 

The ability of the charge distributions on the C^~ 
monomers to rotate creates a picture reminiscent of ro- 
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-25 -20 -15 -10 -5 5 10 15 20 25 
Aw (degrees) 

FIG. 4: Dependence of the potential energy V 1 + V 11 (units 
of kelvin), to which the total potential energy of the crys- 
tal due to C^q-CjTq interactions is proportional, on the de- 
viation angle A-i/'o (units of degrees). The energy scale has 
been shifted for all curves so that the undistorted structure 
(Aipo = 0°) corresponds to zero potential energy. The double 
minimum for KC60 (solid line) shows that an energetically 
more favorable configuration than the Pmnn structure ex- 
ists, featuring rotated electronic densities, while for RbCso 
(dashed line) and CsCgo (dotted line) deviations of the elec- 
tronic distributions on the C<^ monomers from the original 
12 /m structure do not lower the potential energy. 



tating molecules in molecular crystals. It is well known 
that orientational motion of molecular ions in molecular 
crystals influences the translational movements of neigh- 
boring counterions. Here, we have an analogous situation 
and in KCgo one expects average center-of-mass displace- 
ments of the K + ions induced by the angular deviations 
of the charge on the Cg ions. 

The theory of bilinear translation-rotation (TR) cou- 
pling in molecular crystals [E2, 23 is the tool to examine 



the effect of molecular orientations on lattice displace- 
ments of counterions and can be applied here as well. The 
interactions to be considered are the Coulomb attractions 
between the CJ n monomers and the K + ions. Full details 



are given m Appendix |[ The main result, Eq. ( |B15|) , is 
that average displacements (ua(^a)) of the alkali-metal 
ions, labeled by lattice indices t?a = (jiai, ^A2, «A3), are 
found to occur: 



Ui 

(«.M'7 A ))= | | (-1)** 

Here &a is the integer part of tiai + ua3- 
n A i + n A3 = k A + -. 



(6.1) 



(6.2) 



Equation (3.1) shows that there are no displacements 
of the K + ions (and therefore no periodicity doubling) 



along the b axis, and that along the a and c axes dis- 
placements of the K + ions happen in such a way that 
the periodicity doubling scheme of the charge distribu- 
tions on the C^g monomers is respected. The average 
K + displacements — a result of the rotational deviations 
of the charge distributions — constitute therefore a sec- 
ondary doubling mechanism and form a part of the struc- 
tural change in polymerized KCgo- As mentioned in the 
Introduction, displacements of the K + ions were sug- 
gested by Coulon et al. jll| to explain (partly) the struc- 
tural phase transition. The total doubling mechanism, 
now including the K + displacements, is visualized in Fig. 

In Rb- and CsCgo, the charge distributions of the C 60 
monomers do not rotate away from the original 12 /m 
structure. The alkali-metal ions will therefore exhibit 
no average displacements, keeping the 12 /m structure 
intact. 



VII. DISCUSSION AND CONCLUSIONS 

In the preceding sections, we have suggested and ex- 
amined the possibility of having small rotations of the 
valence charge distributions on all C^~ monomers in the 
AC60 alkali-metal fullerides. In KCgo, these electronic 
density distortions lower the potential energy and result 
in average center-of-mass displacements of the K + ions. 
Both the orientational deviations of the charge distribu- 
tions and the translations of the alkali-metals occur in 
such a way that the "new" structure has a doubled peri- 
odicity along the a and c directions in comparison with 
the "old" structure, while the periodicity along the b di- 
rection does not change. These two structural changes 
can therefore account for the experimentally observed 
(a + c,b,a — 2) superstructure in KCgo ■ In RbCgo , the 
potential energy is not lowered by deviations of the ori- 
entations of the charge densities. Hence, the 12 /m struc- 
ture will be preserved. (Since the charge densities do not 
deviate from their 12 /m equilibrium orientations, there 
is no driving force to displace the Rb + ions, which will 
therefore remain at their equilibrium positions.) To sum- 
marize, the model we present forms a possible mechanism 
to explain both the experimentally observed periodicity 
doubling in KCgo and the absence of a similar doubling 
scheme in RbCgo- It establishes a theoretical basis of 
"an appealing hypothesis" discussed in Ref. O where a 
combination of a charge density wave (CDW) with large 
correlated K displacements was suggested as a mecha- 
nism for the superstructure. 

We recall that our model allows rotations of the charge 
distributions associated with the C^~ monomers in the 
lattice, while the C cores remain at fixed positions, i.e., 
the structure formed by the C cores does not deviate from 
the original Pmnn or 12 /m structure. This immediately 
launches the question why the C core network and the 
electronic distribution around it would behave so "inde- 
pendently" (in KC6o). First, we note that the angular 
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FIG. 5: The full doubling mechanism in KC6o: the charge distributions are rotated ±13° away from their original Pmnn 
structure, and the equilibrium positions of the alkali-metal ions are shifted, resulting in a (a + c,b, a — c) superstructure. The 
charge distributions are represented by their contourplot projections of Fig. Q. The radius of the C^" units has been reduced for 
clarity. The alkali-metal ions are represented by filled and empty circles, (a) Projection onto the crystallographic (6, c) plane 
(compare with Fig. H) . (b) Projection onto the crystallographic (a, 5) plane. 



deviation is relatively small: the energy minimum occurs 
at AipQ ~ 13° • Furthermore, one can argue that any de- 
viation angle AipQ different from zero already causes an 
energy lowering, and that it may even be so that the cores 
of the Ceo clusters do "follow" the electronic density dis- 
tortions, thereby causing some restoring forces that pre- 
vent the structure from going as far as the A^o ~ 13° 
configuration but rather causing an equilibrium situation 
at a smaller deviation angle. A smaller equilibrium angle 
can also be the result of restoring forces that act against 
the change in chemical bonding between C^~ monomers, 
since the chemical bonding is affected by rotations of the 
electron distributions on the C^~ monomers. 

Concerning the accompanying metal-insulator transi- 



tion in KCeo jnj, we note that the periodicity doubling 
along the a axis via the mechanism described in Sec. 
[ij will affect the electron transport properties along the 
polymer chains. We point out that the suggested dou- 
bling scheme along the polymerization axis can be seen as 
a CDW, however, not in the usual sense of a charge quan- 
tity modulation, but a modulation of the orientation of 
the charge distribution along the polymer chain, which 
we call an orientational charge density wave (OCDW). 
The charge of every Cgo unit remains the same. It is 
well known that a one-dimensional electron gas, coupled 
to the underlying crystal lattice through electron-phonon 
interactions, is unstable. The Peierls instability leads to 
a CDW accompanied by periodic lattice distortions p3. 
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In the case of a complete charge transfer of one electron to 
each anion, a half-filled band leads to an instability with 
wave vector q = 2fci? = 2{jt/2a), which corresponds to a 
doubling of the unit cell from lattice constant a to 2a in 
real space. The insulating state results from the opening 
of an energy gap that separates the filled lower electron 
band from the empty conduction band. In the present 
case of an OCDW, the modulation of the orientations of 
the charge distributions along the polymer chain and the 
concomitant displacements of the K + ions play the role of 
the lattice distortions, and the metal- insulator transition 
is a consequence of the structural transformation. It is 
not necessarily accompanied by displacements of the C^~ 
monomers along the polymerization direction. Another 
consequence of the rotations of the electron densities on 
the Ceo units is a decrease of the transfer integrals be- 
tween neighboring molecules in a polymer chain, which 
also results in a reduction of the conductivity during the 
phase transition. 

Theoretical work j^, |l^| on the unpolymerized — * poly- 
merized phase transition in the ACgo, A=K, Rb, Cs, com- 
pounds has revealed that the structural difference of the 
polymer phases (space group Pmnn for KCgo, 12 /m for 
Rb- and CsCeo) is due to the electronic quadrupolariz- 
ability of the alkali- metals, and not due to some other 
alkali-mctal-specific parameters such as lattice constants 
or interaction strengths. In our view, the alkali-atom- 
specificity of the structural phase transition (which is 
present for K- but absent for RbCgo) studied here is again 
not due to lattice constants, but is a direct consequence 
of the different space groups of the two compounds stud- 
ied. One can therefore say that it is again the electronic 
quadrupolarizability causing — indirectly — the structural 
difference of K- and RbCgo, since it is responsible for the 
different space groups. 

In conclusion, we have presented a model that (i) ex- 



plains the occurrence of a second structural phase tran- 
sition in KCgo and the absence of such a transition 
in RbCeo, both observed experimentally, and (ii) is a 
starting point for investigations concerning the electronic 
properties of the alkali-metal fullcrides in general and 
the experimentally observed metal-insulator transition in 
KCgo in particular. 
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APPENDIX A: CONDENSATION SCHEME FOR 
THE STRUCTURAL PHASE TRANSITION 

From Fig. |J, it follows that in KCgo, simultaneous de- 
viations of the charge distributions on all C^~ monomers, 
described by the doubling mechanism of Sec. O, lead to 
an energetically more favorable structure. In this ap- 
pendix, we show rigorously that this scheme will indeed 
occur, without making the a priori assumption of alter- 
nating electronic deviations along the a and c axes and 
equal deviati ons along the b axis, expressed mathemati- 
cally by Eq. (£J). 

To each site n, we assign a quantity S(n), which 
takes on the values +1 and —1, corresponding, respec- 
tively, to a deviation angle of the electronic density 
Aip(n) = +Aipo and — A^o- Allowing only the two devi- 
ation angles resulting in minima in the potential energy 
curve (Fig. 0), |A ~0q| = 13.0208°, the interaction energy 
U(n, ft + ft) [Eq. (3.2a)] can be written as 



U(n, n + ft) 



1 + S(n) 



U ++ (n,n + ft) 



1 + S(n + ft) 



S(n) TT „ ,l + S(n + fl) 
— —U + (n,n + fl) — 



■ l±|^^-(n,n + /Q 1 ~ S( " + ^ ) 



S(n) _ 1 -S(n + ft) 
U (n,n + n) , (Al) 



where U ±+ (n,n + ft) is the value of U(ft,n + ft) when 
S(n) = ±1 and S(n+ft) = 1. Analogously, J7 ±_ (n,H+ft) 
is the value of U(n, n+ft) when S(n) = ±1 and S(n+ft) = 
— 1. As in Sec. EO, we consider 14 nearest neighbors. 
The energies [/ ±+ (r?, n + ft ) and C/ ± ~(fi, ft+ft) have been 
calculated using Eqs. (5.1a)-( |5.1c| ) and are listed in Table 
0. Since the corner points (n£ Z 3 ) and the center point 
[n E (Z + i) 3 ] of the orthorhombic cells have a different 
chai n angle (+V , kc 60 anc ^ — V'KCeo' respectively; see Eq. 
( |2.2|) and Fig. |l|), they have to be considered separately. 

The total energy U [Eq. (EO)] is then obtained by sum- 



ming U(ft, H + ft) over the whole lattice. It is convenient 
to write the result in Fourier space. Defining the discrete 
Fourier transform of S(n) by 



S(n) 



_Y^ e ^X(n) S{(f)i (A2a) 



S(q) = —J2e-^^S(n), (A2b) 
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TABLE II: Energies U ++ , U , U + ~ and U , calculated for the 14 nearest neighbors whose relative lattice indices p are 
listed in the first column. A distinction between the corner points (upper part of the table) and the center point (lower part of 
the table) of the orthorhombic cells has to be made. Because of symmetry reasons, only a limited number of different numerical 
values occurs. 
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we obtain 



and 



n jl 

- ^J(9)S( ? 15(- 9 -) + C. 



(A3) 



Here we have split the summation over n into two parts: 
corner points n £ Z 3 and cente r poi nts n E (Z + ^) 3 - 
The summation over n in Eq. (A2b) is understood to 
be a summation over corner points — of which there are 
J\f = N/2 — only. In Eq. (A3), C is an irrelevant constant 



Ja 
J,, 
Jc 
J abc 



J a cos(q x a) + J b cos(q Y b) + J c cos(q z c) 

+ Jal 



+J abc cos(qx^) cos(qy^) cos(q z ^), (A4a) 

. Ja), 

J b ++ -2J b - + J b -, 
J++-2J- + J--, 

^(2J a bc — ^ abc — Jabc )■ 



(A4b) 

(A4c) 
(A4d) 
(A4e) 



The coefficients J+ , . . . , J^ b ~ are related to the poten- 
tial energies U ++ , . . . , U~~ as is indicated in Table |J. 
We now determine the absolute minimum of the func- 
tion J(q) in reciprocal space. Taking into account the 
numerical values of J a = 6 K, J\, = — 14 K, J c = 1 K, 
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and J a bc = 12 K, one finds that the absolute minimum 
lies at q = (qx,QY, 1z) = (f , 0, f ) = qh- The dominance 
of J(qs) leads to a condensation of S(q) at <f = qs- 



(S(q)) = V 



(A5) 



Here rj is the order parameter a mp litude. The conden- 
sation scheme in Fourier space ( |A5| ) corresponds to the 
following real space condensation scheme: 



(S(n))=r ) cos[q B -X(n)}=r 1 (-l) 



(A6) 



APPENDIX B: TRANSLATION- ROTATION 
COUPLING IN KC 60 



where the subsc ript a refers to the alkali-metal ion. Sim- 
ilarly as in Eq. (3.2a), the integration variables r and r*A 
in Eq. (Bl) refer to the local coordinate systems associ- 
ated with the lattice sites H and Ha, respectively, being 
the reason for the appearance of the relative lattice vector 
X(Ha — ft) and the alkali- metal lattice displacement vec- 
tor «a(%)- To be consistent with the earlier convention 



ni 


€ Z ' 






e Z 


< n 2 


e Z 


H 




e Z 


> "3 


e Z t 




k "3 


e Z 



(B2) 



one must have for the alkali- metal lattice indices ft a- 



In this appendix, we examine the coupling between the 
orientational deviations of the charge distributions on the 
C^~q monomers — lowering the crystal's potential energy in 
KCgo (see Fig. ||) — and displacements of the K + ions. We 
use concepts of the theory of bilinear translation-rotation 
(TR) coupling in molecular crystals 22, ^3), which is 
generally used to determine the influence of molecular 
orientations on lattice displacements of counterions. 

The starting point is the potential energy of a C^~ 
monomer and a K + ion: 



U (H, n A ) 



1 



dr I drA 



4ire 

p(f; tp(H) + Aip(H))p A (r A ) 
\r- r A - X(H A - ft) - 4(4) 



(Bl) 



UAl S Z 
UA2 € Z 
n A3 e Z - 




(B3) 



We treat the K + ion as a point charge and write for 
its charge distribution: 



Pa(ta) = eS(f A )- 



(B4) 



Working with the charge distribution p(r; ip(n) + 
Aip(H)) = p(r,9,(j)- ijj(ft) ~ Aijj(ft)) of Sec |y| for the 
C^q monomer and taking into account Eq. QB4| ), we get 
the following expression for U(H, Ha)'- 



U(H,H A ) 
d(9,4>;H,H A ) 



F I sui 
'o 



2tt 



,p a {e,4>-i>(n)- A^(rt)) 



d(6, <f>] ft, Ha) 

{[Rcos9 — (uai — ni)a — uai("a)] 2 + [i? sin cos </!> — (nA2 — n 2 )b — ua2(Ha)Y 

1 /2 

-I- [R sin 6 sin <j>- (n A3 - n 3 )c - ua3("a)] 2 } 



(B5a) 



(B5b) 



where the constant F is given by Eq. ( 5.11 ). 

We consider small center-of-mass displacements of the 
K + ions and expand U (ft, Ha) in terms of the compo- 
nents of 4 (ft a), retaining only the zeroth- and first-order 
terms: 



±1, corresponding to Aip(H) — ±|A?/>o|, with |A^o| = 
13.0208°. The interaction energy U(H,Ha) can then be 
written as 



U(ft,ft A ) = U(n,n A )\ aA{iiA)= o 



dU(H, Ha) 



y 

f-f du A i(n A ) 



uai(Ha)- (B6) 



«a("a)=0 



As in Appendix [A|, we introduce the quantity S(H) = 



U(H,H A ) = 1 -±^U+(H,Ha) + 1 -^±U-(H,Ha), 

(B7) 

where U^(H,Ha) is the value of U (ft, ft a) when S(n) = 
±1. The expansion (B6) becomes 
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1 ~\~ S(fi) I 

U(n,n A ) = < V + {n,n A ) + ^v+ (ri,ri A )u Al (n A ) 

{ i=l 

1 — S(n) _ . _,_ _ . , 

H ^ p (n,n A ) + > v,, (n,n A )u Al ( 

r 



HA J 



(B8a) 



with 



V±{n,n A ) = U^ft^^, (B8b) 
dU ± (n, n A ) 



v i (n,n A ) 



du Ai (n A ) 



«a("a)=0 



c) 



actions to the potential energy is obtained by summing 
U (ft, ri A ) over the whole crystal lattice: 



The contribution Uc &0 -A of all electrostatic C 60 -K + inter- where 



W C60 -A=W R +W TR , 



(B9a) 



U K = 



U 



EE 

n (3a 

-EE 



n + /i A )H K (n,n + ^ A )^, 



■ 2J «<"(ra, u + fl A )u Ai (n + p A ) + 



(B9b) 

^v^(ri 7 n + fl A )u Ai (ri + fl A ) \ . (B9c) 



For a given C^~ monomer, we limit ourselves to the six where m A is the mass of the K + ion, and using in addition 
nearest alkali- metal neighbors. The values of y ± (n, n+ t he d iscrete Fou rier transform of S(n), defined by Eqs. 
u A ) and v^(n,H + j2 A ), i — 1,2,3 are listed in Table ( A2a ) and ( A2h| ), we get for the TR term of the potential 
H. As in Appendix |A|, a distinction has to be made energy 
between the corner points (n € 2 3 ) and the center point 
[n E (Z + i) 3 ] of the orthorhombic cells. 

Introducing the discrete Fourier transforms of the 
alkali-metal ion displacements, 



U TR = 



u A (n A ) 



1 



im A 



J2S(-q)v(q)-u A (q), (Blla) 
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(Bllb) 



and In the theory of bilinear TR coupling in orientationally 

disordered crystals |23|, , it is shown that the minimal 
potential energy for a given orientational configuration 
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TABLE III: Energies V* and derivatives of energies v^, i = 1,2,3, calculated for the six nearest alkali-metal neighbors with 
relative lattice indices Pa- A distinction between the corner points (upper part of the table) and the center point (lower part 
of the table) of the orthorhombic cells has to be made. Only a limited number of different numerical values occurs. 
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K A^ 1 


+^2 + c 




+< 



{S(q)} is obtained when 



U A (q) = -M-\q)v(-q)S(q), 



(B12) 



where M{q) is the dynamical matrix of the orthorhom- 
bic crystal in absence of TR coupling. Using for S(q) 
the minimal po tent ial energy condensation scheme of Ap- 
pendix Eq. (A5), we find for the displacements in re- 
ciprocal space: 



/ -42? 



(u A (0)) =-M-\q = q B )i 



lab 



I WNSqjB- 

(B13) 

The numerical values of 2?i a h and 2?3 C follow from Table 
iTfl: V lab = -49.7 K A" 1 and 2? 3c = 222.9 K A -1 . Since 



qs = (— ,0, -), the scheme (B13) becomes in real space 
1 



(u A (n A )) 



-42? lai , 
M- l {q = q B ) | | 'M 
-22? 3c 



Xe -^("Ai+nA 3 )_ ( B14 ) 

The alkali-metal ions are located on a sublattice obeying 
?iai+«A3 = ^A+l with k A 6 Z [seeEq. (B3)]. Therefore, 
one has e -™(™^i+«A 3 ) = _,;(_i)feA ; yielding 

/ -42? lQb 

I r?(-l) fcA 



1 



(u A (n A )) = -=M (q = q B ) 




\ -22? 3c 



(B15) 
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